Klein tunneling of light in fiber Bragg gratings 

Stefano Longhi 
Dipartimento di Fisica, Politecnico di Milano, 
Piazza L. da Vinci 32, 1-20133 Milano, Italy 

Abstract 

A photonic analogue of Klein tunneling (KT), i.e. of the exotic property of relativistic electrons to 
pass a large repulsive and sharp potential step, is proposed for pulse propagation in a nonuniform 
fiber Bragg grating with an embedded chirped region. KT can be simply observed as the opening of 
a transmission window inside the grating stop band, provided that the impressed chirp is realized 
over a length of the order of the analogue of the Compton wavelength. 
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I. INTRODUCTION 



A remarkable prediction of the Dirac equation is that a below-barrier electron can pass 
a large repulsive and sharp potential step without the exponential damping expected for a 
non-relativistic particle. Such a transparency effect, originally predicted by Klein jlj and 
referred to as Klein tunneling (KT), arises from the existence of negative-energy solutions 
of the Dirac equation and requires a potential step height AV of the order of twice the rest 
energy mc 2 of the electron [2j. Relativistic tunneling across a smooth potential step, which 
describes the more physical situation of a constant electric field E in a finite region of space 
of length I, was subsequently studied by Sauter j^]. Sauter showed that to observe barrier 
transparency the potential increase AV ~ eEl should occur over a distance I of the order or 
smaller than the Compton wavelength Ac = h/(mc), the transmission probability rapidly 
decaying toward zero for a smoother potential increase 2j-|4|. The required field corresponds 
to the critical field for e + e~ pair production in vacuum, and its value is extremely strong 
making the observation of relativistic KT for electrons very challenging. Therefore, growing 



efforts have been devoted to find experimentally accessible systems to investigate analogs of 



relativistic KT 



a. 



Recently, great interest has suscitated the proposal [6] and first exper- 
imental evidences of KT for non-relativistic electrons in graphene, which behave like 
massless Dirac fermions. On the other hand, optics has offered on many occasions a test 
bed to investigate the dynamical aspects embodied in a wide variety of coherent quantum 
phenomena (see, for instance, js] and references therein). In optics, several proposals of KT 
analogs have been suggested as well, including light propagation in deformed honeycomb 



photonic lattices 



10( whose band structure is similar to the one of graphene ll|, |12J], light 



13], spatial light 



refraction at the interface between positive-index and negative-index media 
propagation in binary waveguide arrays [l^] , and stationary light pulses in an atomic ensem- 



ble with electromagnetically induced transparency [15J. The experimental implementations 
of such schemes, however, might be a nontrivial matter, and an experimental observation of 
KT for photons is still lacking. On the other hand, multilayer and Bragg dielectric struc- 
tures, such as fiber Bragg gratings (FBGs), are rather simple photonic devices with flexible 
design that have been successfully demonstrated to provide an accessible laboratory tool to 



investigate photonic analogues of non-relativistic tunneling phenomena 



lj|-18|. Here it is 



shown that an optical analogue of KT can be achieved in a nonuniform FBG composed by 
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two periodic sections linked by a chirped section which mimics an external potential step 
in the Dirac equation. Such a FBG-based system might be considered the simplest system 
proposed so far in order to observe Klein tunneling in any optical system. 



II. QUANTUM-OPTICAL ANALOGY 

The starting point of our analysis if provided by a standard model of light propagation 
in a FBG with a longitudinal refractive index n(z') = n + An m(z') cos[27r;z'/A + 2<p(z')], 
where n is the effective mode index in absence of the grating, An <C n Q is the peak index 
change of the grating, A is the nominal period of the grating defining the reference frequency 
ujb = vrc/(An ) of Bragg scattering, c is the speed of light in vacuum, and m(z'), 2cj)(z') 
describe the slow variation, as compared to the scale of A, of normalized amplitude and 
phase, respectively, of the index modulation. Note that the local spatial frequency of the 
grating is k(z') = 2tt / A + 2(d<f)/dz'), so that the local chirp rate is C = dk/dz' = 2(d 2 <f)/dz' 2 ). 
The periodic index modulation leads to Bragg scattering between two counterpropagating 
waves at frequencies close to u>b- By letting E(z', t) = (fi(z', t) exp[—icoBt + iksz' + i(f>(z')] + 
ipi{z' , t) exp[— iust — iksz' — i(j){z')\ + c.c. for the electric field in the fiber, where fee = 7r/A, 
;he envelopes tp\ and ^ 2 of counterpropagating waves satisfy the coupled-mode equations 



where k(z') = [kBm(z')An]/(2no) and v g ~ c/uq is the group velocity at the Bragg frequency. 
The analogy between pulse propagation in the FBG and the Dirac equation in presence of 
an electrostatic field is at best captured by introducing the dimensionless variables z = z'/Z 
and r = t/T , with characteristic spatial and time scales Z = 2no/(&BAn) and T = Z/v g , 
and the new envelopes vpi^iz') = [(pi(z') =F ¥2(2')]/ V%- In this way, Eqs.(l-2) can be cast in 
the Dirac form 




i [<9 Z / + {l/v g )d t \ ipi = (d(f)/dz')(pi - K(J)<p 2 
i[-d X ' + (l/v g )dt](p2 = {d<p/dz')(p 2 - k(z')lpi 



(1) 
(2) 



id T ip = —iaid z ip + ma 3 i/; + V(z)ip 



(3) 
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for the spinor wave function ip = (ipi, ip2) T , where V(z) = (d<p/dz) and 01,3 are the Pauli 
matrices, defined by 

'\ 
-1 




(4) 
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FIG. 1: Energy diagrams of the Dirac equation (3) for (a) a sharp, and (b) a smooth potential step 
V(z) of height AV. The shaded regions are the forbidden energies that separate the electron and 
positron states, the dotted horizontal line is the energy of the incoming wave packet, and the 
dashed curve is the shape of the potential step V(z). (c) Schematic of the grating structure that 
realizes the optical analogue of relativistic tunneling across a potential step. The grating comprises 
five sections, denoted by roman numbers (from I to V), that are defined by the amplitude m(z) 
and phase gradient V(z) = (d(j)/dz) grating profiles. 

In its present form, Eq.(3) is formally analogous to the one-dimensional Dirac equation 
with h — c — 1 in presence of an external electrostatic potential V(z), m playing the role of a 
dimensionless (and generally space-dependent) rest mass (see, for instance, |2|,|4|). As is well- 
known, a nonvanishing mass m is responsible for the existence of a forbidden energy region, 
which separates the positive- and negative-energy branches of the massive Dirac equation. 
The optical analogue of the forbidden energy region is precisely the photonic stop band of 
the periodic grating. As the refractive index modulation of the grating, i.e. the mass term 
m in the Dirac equation (3), is decreased, the stop band region shrinks and the limit of a 



massless Dirac equation (similar to the one describing the dynamics of electrons in graphene 
near a Dirac point) is attained. The additional external potential V in Eq.(3), related to the 
chirp of the grating according to V(z) = (d(p/dz), changes the local position of the forbidden 
energy region. Therefore, pulse propagation in a FBG with a suitably designed chirp profile 
can be used to mimic the relativistic tunneling of a wave packet in a potential step V(z). 
It should be noted that, as compared to other photonic analogues of KT recently proposed 
in Refs. jlo|, Q and based on spatial light propagation in periodic photonic structures, the 
phenomenon of KT occurring in FBGs and discussed in the following section involves the 
temporal (rather than the spatial) light dynamics and can be therefore simply investigated 
in the frequency domain by spectrally-resolved transmission measurements. 

III. KLEIN TUNNELING 

To realize the analogue of KT, let us first assume that the optical pulse propagates in 
a region of the grating where m(z) is uniform and equal to one, and let us assume a chirp 
profile that mimics a step potential with an increase from V = to V = AV which occurs 
over a length I (see Fig.l). Since for the Dirac equation (3) written in dimensionless units 
the Compton length is Ac = 1 and the rest energy is mc 2 = 1, according to Sauter's analysis 
KT is expected to be observable for I smaller than ~ 1 and for a potential height AV 



larger than 2 



The process of KT and tunneling inhibition for a smooth potential 



step can be simply explained by a graphical analysis of the space-energy diagrams (z, Q) of 
the one-dimensional Dirac equation 2j, which are shown in Figs. 1(a) and (b) for a sharp 
and for a smooth potential step, respectively. For the sake of clearness, in the figures the 
potential V(z) has been chosen to yield a nonvanishing and constant chirp rate over a length 
/; different forms for the potential step, such as the profile V(z) = (AV/2)[1 + tanh(z//)] 
considered in the seminal work by Sauter [3j, can be assumed as well without changing 
the main results. The space-energy diagrams of Figs. 1(a) and (b) schematically show the 
behavior of the energy spectrum of Eq.(3) versus z, which is composed by two branches 
-the electron and positron energy branches of the Dirac equation- separated by a gap of 
width 2m(z) and centered along the curve Q = V(z). The gap regions are visualized in 
the diagrams by the shaded areas. A wave packet (optical pulse) in the electron branch 
with an initial mean energy Q (1 < Q < AV — 1) coming from z — > — oo tunnels into the 
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FIG. 2: (a-c) Pulse propagation in a FBG with a chirp profile V(z) = (AV/2)[1 + tanh(z/Z)] for 
= 6 and for: (a) I = 0.1, (b) I = 1.5, and (c) Z = 5. (d) Profiles of V = (d<p/dz) (upper 
plot) and of grating amplitude m(z) (apodization profile, lower plot). The potential V is shown 
for I = 5, corresponding to the simulation of Fig.2(c). 



z > region after crossing a forbidden energy region, indicated by the bold segment AB 
in Fig. 1(b), which vanishes for a sharp potential step [I = 0, see Fig. 1(a)]. According to 
Sauter's analysis , the tunneling probability is appreciable provided that I is smaller 
than ~ 1. In the FBG context, the energy diagrams of Fig.l are equivalent to the band- 
reflection diagrams introduced by Poladian for a graphical analysis of nonuniform gratings 
[20I . where the energy represents the frequency detuning of the incoming wave from the 
Bragg frequenc y cu r. The Sauter's condition I <~ 1 for KT can be derived following the 
analysis of Ref. 20| by computation of the transmittance of the effective grating associated 
to the evanescent region AB shown in Fig. 1(b) (see Sec.V.A of Ref. [20I] ) . In the previous 
discussion, we assumed m(z) = 1, however for a grating with finite spatial extent one has 
m(z) — > as z — > ±00. To inject and to eject the optical pulse into the m(z) = 1 grating 
region around z = 0, an input and an output apodization sections can be introduced, which 
adiabatically convert the input and output wave packets from the m(z) = regions into 
the m(z) = 1 grating region (see Fig.l). Therefore, the general structure of the FBG that 
realizes a photonic analogue of relativistic tunneling across a potential step consists of five 
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sections, as shown in Fig. 1(c): two boundary apodization sections (regions I and V), and 
two uniform sections (regions II and IV) separated by a central chirped section of length 
~ I (region III). In Figs. 2(a-c) typical examples of pulse tunneling across the potential 
step V(z) = (AV/2)[1 + tanh(z/Z)] are presented, showing KT for a sharp potential step 
[Fig. 2(a)] and inhibition of tunneling as the step gets smooth [Fig.2(b) and (c)]. The figures 
depict the temporal evolution of \ipi\ 2 + l^l 2 = Iv^il 2 + Iv^l 2 -which is proportional to the 
field intensity averaged in time over a few optical cycles and in space over a few wavelengths- 
as obtained by numerical analysis of Eqs.(l) and (2) for a grating length of z = 160 with 
a quarter-cosine apodization profile [see Fig.2(d)], AV^ = 6, and for a few values of I. A 
forward-propagating Gaussian pulse ipi of mean energy Q = 2 coming from z — > — oo and 
of duration (FWHM in intensity) t p = 5 has been assumed as an initial condition. For 
typical parameter values no = 1.45, An = 3.3 x 10~ 4 and Xb = 2^oA = 1560 nm, which 
apply to FBGs used in optical communications, the spatial and temporal scales in Fig.2 are 
Z ~ 1.5 mm and T ~ 7.3 ps, respectively. Hence, in physical units the grating length is 
L ~ 24 cm, whereas the optical analogue of the Compton length is Ac = Z ~ 1.5 mm. 
Such nonuniform FBG structures should be realizable with current FBG technology based 



on UV continuous laser writing 



21] . It should be finally noticed that, as in an experiment 



the tracing of pulse evolution in the grating (Fig.2) can be a nontrivial task, the signatures 
of KT can be simply obtained from standard spectral transmission measurements of the 
grating. In fact, for a given value of AV > 2 and according to the band diagram of Fig.l(a), 
in the KT regime a transmission window at 1 < Q < AV — 1 [the segment AB in Fig. 1(a)], 
embedded into the two gaps A^ — 1 < Q < AV + 1 and — 1 < Q < 1 [the segments BC and 
AD in Fig. 1(a)] should be observed in the transmission spectrum, the suppression of KT for 
a smooth potential corresponding to the lowering of such a transmission window. This is 
clearly shown in Fig.3, where the spectral transmittance of the FBGs corresponding to the 
simulations of Figs. 2(a), (b) and (c) are depicted. Note that, as the length I of the chirped 
region is increased [from Fig.3(a) to 3(c)], the transmission window embedded in the two 
adjacent gaps disappears, which is the signature of KT inhibition. 
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FIG. 3: Numerically computed spectral transmittance of FBGs used in numerical simulations of 
Fig. 2 for no = 1.45, An = 3.3 x 10~ 4 , Xb = 1560 nm, corresponding to a grating length L ~ 24 
cm, and for increasing values of I: (a) I = 0.1, (b) I = 1.5, and (c) / = 5. 

IV. CONCLUSIONS 

In conclusion, a photonic analogue of Klein tunneling based on pulse propagation in 
nonuniform fiber Bragg gratings has been proposed. As compared to other photonic ana- 



logues of KT recently proposed in Refs. 



10 



141 ] and based on spatial light propagation in 



periodic photonic structures, the phenomenon of KT in FBGs suggested in this work can 
be simply observed in the frequency domain as the opening of a transmission window inside 
the grating stop band, provided that the impressed chirp is realized over a length of the 
order of the analogue of the Compton wavelength. Such a FBG-based system might be thus 
considered to be the simplest optical analogue proposed so far to observe KT. 
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